Friedel oscillations of Kondo impurities: A comparison. 
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Abstract 

Recently Affleck et al. derived the existence of Friedel oscillations in the presence 
of a Kondo impurity. They supported their analytic derivation by numerical calcu- 
lations using Wilson's renormalization approach (NRG). In this paper the size of the 
Friedel oscillations is calculated with the FAIR method (Friedel Artificially Inserted 
Resonance) which has been recently developed. The results of NRG and FAIR are 
compared. The development of Friedel oscillations with a phase shift of 7r/2 outside of 
the Kondo radius is confirmed. 



The properties of magnetic impurities in a metal is a fascinating problem which was first 
studied by Friedel [1] and Anderson [2] . The disappearance of the magnetic moment at low 
temperatures, the Kondo effect, is one of the most intensively studied problems in solid state 
physics [3], [4], [5], [6], [Tj, [8], [9], [lOj, [ll], [12], [l3], [M], [15], [l6], [I7], [H]. In the last 
decade the Kondo effect has experienced a renaissance. There is a growing interest in this 
field [I9], extending from magnetic atoms on the surface of corrals [20] to carbon nanotubes 
[21], quantum dots [22], [23], [21], [25], [26], [27], [28] and nanostructures [29]. There are still 
many open questions, particularly the real-space form of the wave function and the resulting 
charge density and polarization. 

Recently Affleck, Borda and Saleur [30] (ABS) investigated the formation of Friedel 
oscillations in the vicinity of a Kondo impurity. Their result has the form 
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where D is the dimension of the system, the coefficients Cd have the values Ci = 1/ (27r), 
6*2 = 1/ (27r^) and C3 = 1/ (47r^) , and = hvp/kBTx is the Kondo length. The function 
F{r/rK) is a universal function which approaches the values +1 for r/r^ << 1 and —1 for 
r/rx » 1. (I skipped the phase shift 6p due to potential scattering). 
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Besides the analytic derivation of ([T]) ABS performed also numerical calculations using 
Wilson's NRG approach [TT]. In the NRG technique one uses Wilson states with a loga- 
rithmic discretization of the conduction band (which extends from (—1, +1) (in units of the 
Fermi energy) and has a constant density of states.) 

ABS point out that Wilson energy states are not well suited for the calculation of such 
fine effects as Friedel oscillations. The reason is that the Wilson states average over a large 
number of original k-states. Therefore they use a modified method which was introduced ear- 
lier by Borda (for details see [2n])- With this approach they obtain a numerical confirmation 
of the universal function F{r/^x)- 

The author has introduced during the past few years a new numerical approach to the 
Kondo and the Friedel-Anderson impurity, the FAIR-method (Friedel Artificially Inserted 
Resonance). It is based on the fact the n-electron ground state of the Friedel Hamiltonian 
(consisting of an electron band and a d-resonance) can be exactly expressed as the sum of 
two Slater states 

n— 1 rt— 1 

= Aal n 4% + Bd^ n 4'^o (2) 

i=l i=l 

where aj is an artificial Friedel resonance state which determines uniquely the full orthonor- 
mal basis l*^!!- extension of this ground state to the Friedel-Anderson and Kondo 
impurity yields good numerical results. Recently this method was applied to calculate the 
Kondo polarization cloud for those impurities [31]. Therefore the author could not resist 
using the FAIR method to calculate the Friedel oscillations of a Kondo impurity. As in the 
paper by ABS I am treating the one-dimensional case. (The two- and three-dimensional 
treatments are essentially identical). Also the calculation is restricted to zero impurity scat- 
tering (the same as in ABS's calculations). 

The construction of the Wilson states is essentially the same in NRG and FAIR. One 
starts from an electron band with linear dispersion as shown in Fig.l. If one follows Wilson 
and measures the energy in units of the Fermi energy and the wave vector k, = k/kp in units 
of kp then the dispersion is just 

e^ = K-1 

This band contains a macroscopic number of states (f^. (~ 10^^). Wilson subdivided 
the positive and negative half of the energy band by a factor A. For example for A = 2 
the negative band is divided into energy cells Co = (—1 : — |) , = (— | : —\) ,■■ = 
(Cz/ : Cj^+i), where = —2'". In each energy cell (t^ one has Z^, states ip^. with the energy 
{k — 1). Assuming that each state ip^ [f) has the same interaction with the impurity then 
the states in the cell (t^, are combined to a new state ipi, (r) 

(r) = (r) 

V 

This state represents the full interaction of all states (r) in <ti, with the impurity. (Its 
interaction is enhanced by \fZl\ The sub-Hilbert space of Cj, contains [Zy — 1) additional 
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states 

with / and /i running from to [Z^ — 1) (Here ip^ is identical to ipi,). The states with / > 
have no interaction with the impurity, but they have finite matrix elements |ifo I "V^^/ ) 
with the band energy Hamiltonian Hq. These matrix elements are neglected in the NRG 
and the FAIR calculation. This is strictly speaking only correct when all the original states 
ifi^ (r) in a given energy cell 'ty have the same energy. In Fig.l the energy - wave vector 
diagram is shown. The dashed line gives the linear dispersion relation. The darkly shaded 
square marks the energy cell Cq- The thick zig-zag line gives the dispersion for which the 
Wilson NRG would be exact. 




Fig.l: The effective dispersion relation between energy e and 
wave number k in Wilson's Hamiltonian. 

In the one-dimensional case the wave functions ip^ are cosine and sine functions. (Only 
the cosine states interact with the impurity). Since the energy C, and the wave numbers k 
are given in normalized units oi Sp and kp it is natural to measure distances also in reduced 
units. Throughout this paper lengths are measured in units of Xf/2, for example 
the distance from the impurity is denoted as ^ = r / {Xp/2). 

For the wave functions I use (p^ (0 = v^2/A cos (ttk^) in the range from —A/2 to +A/2. 
(Here A is given by the size L of the one-dimensional box, A = L/ {Xp/2) . The number of 
states ipf, (0 ill the energy cell (ty is {(u+i — Cu) / (2/A) 
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(0 = ... ' ..^sin ( vre^^^±^l cos fvre f 1 + ^^^±i±^^) ) (3) 



This yields for the wave function of the Wilson state ipu (^) 

2 1 

The exchange interaction is given by Hgd 

Hsd = VaJT.a,p'^i (0) (0) ■ S (4) 

where Va is the atomic volume and J is the exchange interaction 
The FAIR method yields for the Kondo ground state 
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As^dal^^d\ + Ad,sd]al__^ |0) + Ad,sal_^d\ + As4d]al^^ |0) 

n— 1 71—1 n—1 71—1 

io) = n «It n «]-i<f io) = n «Lt n 4+1^^ 

i=l j=l i=l j=l 

The two FAIR states aQ_^_ and aj_ are obtained by optimizing the ground-state energy 
(see for example [32], [31]). They have the form 

where cj, are the creation operators of the Wilson states ip^,. (These two states are maximally 
coupled to the impurity. I expect that they are closely related to Wilson's state /o). The 

FAIR states al_^ determine uniquely the two bases The total charge density of the 

ground state can be calculated from the ground state with the help of the individual 
wave functions of the Wilson states. 

To calculate the charge density of the Kondo ground state one has to calculate the charge 
densities of all states al_^ for < z < A^/2 and add them according to the occupation of these 
states. This yields the main contribution. But there are in addition two interference terms 
which yield rather small contributions. The total charge density is then calculated at average 
distances from the impurity of 2, 4, 8, ..2', ...2"^^. This is done in intervals of two wave lengths 
(A^ = 2 or Ar = Xp) where a dense trace of the charge density as a function of ^ is 
calculated. For each trace the position of the charge minima and maxima is determined as 
well as the amplitude. 

In Fig. 2 the numerical result for F (^/^k) is plotted. It is the amplitude of the first term in 
equ. i.e. Cd cos {2kFr - Dfj jr^ (which in one dimensions is equal to sin (27r^) / (27r^)). 
The abscissa is the logarithm (base 2) of the average distance ^ from the impurities. (The 
change in sign of F {C,/C,k) corresponds to a jump in the phase by vr). 

The top curve with open circles is for = 60 Wilson states with a ratio A = 2 between 
neighboring energies ^i^. It confirms the statement by ABS that these Wilson states are not 
well enough suited for the calculation of small charge fluctuations on a length scale of the 
Fermi wave length. 
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To improve the calculation I tried a different path than ABS. In the FAIR method 
the two essential states are the FAIR states aj+and aj_. When they are known one can 
easily construct the full orthonormal basis for each FAIR state (with the condition that 

(^al^o\Ho\aj^o^ — for i,j > 0). This is a trivial calculation with negligible computer 
time. Each aj is composed of Wilson states 

EN-l V t 

where cj, is the creation operator of the Wilson state The square IoiqI gives the occupa- 
tion or contribution of the state cj, or to aj. Since the state is constructed from the orig- 
inal states the ratio of |q;q |^ divided by width of the energy cell, i.e. p,^ = |^ / (Ci/+i — C,v) 
gives the relative contribution of the original states to the FAIR state aj. This state den- 
sity corresponds to a function p (C) with the condition p (C) = 'p^ for Ci/ < C < C+i- The 
integral of p (C) over the energy is normalized 

For large N the state density p {() shows a smooth dependence on the energy. (This 
is shown in Fig. 4 for the example of a Friedel resonance, see below). Therefore one can 
extrapolate p {() from a finite representation N and construct the FAIR state for twice the 
number of states. The most natural way to double the number of Wilson states is to set a 
new energy ratio Cl/Cl+i = A' = \/A. In this way the FAIR states and their full bases are 
obtained for N = 120 (with a new A = ^/2). It turns out that the new solution for 
is better than the old one (for N = 60) because its ground-state energy is lowered. The 
second curve in Fig.2 (open triangles) shows the numerical result for the function F {^/^k)- 
It clearly improves the result although it does not yet reach the value —1 for large ^. A 
second doubling of N yields = 240 and A = \/2. The ratio: width of an energy cell 
divided by the average energy is given by 2 (A — 1) / (A + 1) . For the usual Wilson states 
with A = 2 this ratio is 2/3. After two additional splitting of the states this ratio becomes 
about 0.173. Now the energy states are much closer together. For N = 240 the numerical 
result for F {^/^k) is shown in the bottom curve of Fig.2. It finally shows the two hmiting 
values of at short distances and —1 at large distances in agreement with the theoretical 
prediction. 
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Fig.2: The numerical result for the function F {^/^k) in equ. ([T]) 
for different ratios of neighboring energies Ci/- For the number of Wilson 
states = 60 the ratio is A = 2, for = 120 it is A = ^2 and for N = 240 
it is A = 

Since the function F {^/ C,k) is a universal function a comparison of Kondo impurities with 
different Kondo energies is appropriate. ABS performed a large number of calculations. I 
want to show that the FAIR method yields the same universality. In the past we investigated 
Kondo impurities with different exchange interactions and therefore different Kondo energies 
|32j . Here I consider the two examples with J = 0.1 and J = 0.08. One obtains the Kondo 
lengths from the corresponding the Kondo energies Ek = 2.35 x 10^^ and 1.37 x 10~®. 

In Fig.2 we plotted the amplitude of the first term in equ.([T]) to extract the function 
F {^/^k)- Now we determine the actual amplitude A of Friedel oscillation which 

is given by 

27r4 27r4 

In Fig. 3 this amplitude A {^/^k) = [1 — -F {^/C,k)] is plotted for the two different Kondo 
energies as a function of log2 / C,k) ■ Since the abscissa has a logarithmic scale the curves 
are just shifted by log2 {^k) ■ Although the Kondo energies differ by roughly a factor of five 
the two curves are essentially identical. Therefore the amplitude of the Friedel oscillation is 
universal (for the symmetric case and a linear dispersion relation). The numerical calculation 
yields a phase shift of 7r/2 for distances larger than 2 Fermi wave lengths. For most of the 
region {Q < ^/^k ^ 16) the relative deviation is less than 10~^. 

ABS give two asymptotic function for ^ << and ^ >> C,k- These functions which are 
shown as dotted curves in Fig.3, have been shifted by 0.5 to the left (see upper scale). This 
means that I use a slightly different definition of the Kondo energy than ABS, which yields 
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a different Kondo length. The ratio between our Kondo energies is 2°'^ ^ 1.4. This is not 
surprising because ABS use Wilson's definition which requires the NRG calculation whereas 
we define the Kondo energy as the energy difference between triplet and singlet state [52] . 
ABS obtained from their numerical calculation that the function F{^/^k) crosses zero at 
{^/^k) ~ 0.12 ± 0.02. In Fig.3 this point lies at 2-2-9 = 0.13. 
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Fig.3: The amplitude of the Friedel oscillation [1 + F {^/C,k)] is 
plotted versus ^/^k for two different Kondo energies. 

It is worthwhile mentioning that in this calculation the Kondo ground state is calculated 
once. From this ground state one obtains the Friedel oscillation for all distances ^. ABS had 
to perform a new NRG calculation for each individual distance ^ (which actually represent 
states with slightly different ground-state energies). 

One essential step in calculating the Friedel charge oscillations was the repeated doubling 
of the number of Wilson states. This was possible because the compositions of the FAIR 
states aQ_,_ and al_ are essentially given by the state density functions p+ {() and p_ {(). This 
shall be demonstrated for the FAIR state of a Friedel (resonance) impurity. 

We take the Wilson electron band ranging from (—1 : +1) with constant density of states 
(equal to 1/2) and a d resonance at the Fermi energy {Ed = 0) with a s-d matrix element 
\Vsdf' = 0.1. The ground state is given in ([2]) with the FAIR state aj. 

In Fig. 4 the values p,^ = {a^f / {(u+i — Cu) are plotted at the discrete energies (C+i + Cu) /2. 
The full circles are the values for N = 24, A = 2. Only a small fraction of the whole energy 
band is included in the figure to show an optimal section. The full triangles show the values 
for = 48, A = while the stars represent the values of p,y for = 96, A = v^- For each 
A^ the optimal ground state has been derived. One recognizes that the points follow indeed 
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a common function p(C)- This justifies the doubhng of the Wilson states since the new pi, 
values can be obtained by inter- and extrapolation. This yields the composition of the new 
FAIR state after the doubhng of the number of Wilson states N. 
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As an after thought of this consideration I calculated the Friedel oscillation for Friedel 
impurities with the d-resonance energy at the Fermi level, Ed = 0, and very small resonance 
width A = 71 \ Vsd\^ p where p = 1/2 represents the density of states for the Wilson spectrum. 
By using rather small values for |Kd|^= 0.5 x 10~^, 1.0 x 10~^ and 2.0 x 10~^ one can evaluate 
the amplitude of the Friedel oscillations as a function of ^. 

The numerical calculation yields a charge oscillation Appr = ^Fr (0 cos (27r^ — 5pr) / (Stt,^). 
The phase shift dpr in the range 3 < log2.^ < 19 is equal to 7r/2 within 1% accuracy. For 
shorter distances shorter than the Friedel coherence length {^p^. « 2/ iTry^rfl^) the amplitude 
is strongly suppressed. The amplitude is also universal since the three curves in Fig. 5 would 
perfectly coincide if plotted as a function of (^/^Fr-)- This behavior is quite analogous to 
the Friedel oscillations of a Kondo impurity. However, the functional dependence and the 
long-distance amplitude differ. 
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Fig. 5: The amplitude of the Friedel oscillation for Priedel 
impurities with s-d hopping matrix elements of 
\Vsd\^ = 2 X 10-^ 1 X 10"^ and .5 x 10"^ The Friedel 
oscillations develop only for distances larger than the 
Friedel length ^pr- 

To conclude, the goal of this paper has been to reproduce ABS result for the Friedel 
oscillations of a Kondo impurity with the FAIR method. This is a rather delicate problem 
because electron densities have to be calculated on a spatial scale much smaller than the 
Fermi wave length Xp over distances of up to IQ^Xp with a relative accuracy of A^/r, i.e. 
10~^ or better. The numerical results of the FAIR method have about the same quality as 
the NRG method. 
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